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5.  Report  of  work  done 

Our  work  involves  identifying  light  weight  adsorbates  other  than  metal  atoms  that  are 
used  to  functionalize  CNT  to  enhance  its  storage  capacity.  So  we  have  employed  light 
weight,  hydrogen  rich  metal  and  non-metal  hydrides.  The  total  energy  calculations  are 
carried  out  using  density  functional  theory  implemented  in  VASP.  Ultra  soft 
pseudopotentials  are  used. 

(i)  Study  of  hydrogen  storage  in  single  walled  carbon  nanotube  functionalized  with 

NH3 


Ammonia  or  Azane  exists  in  gaseous  form.  It  has  a  lone  electron  pair  and  so  we  can 
say  it  is  electron  rich.  Thus  our  study  enables  to  understand  the  hydrogen  storage 
capability  of  a  electron  rich  and  neutral  molecule,  NH3.  The  simulation  cell  of 
dimensions  24  A  x  24  A  x  4.98  A  with  two  unit  cells  of  SWCNT  (5,  5)  is  considered. 
The  radius  of  the  (5,  5)  CNT  is  3.44  A.  The  investigation  is  done  for  two  different 
coverages  namely  half  and  full  coverages.  Half  coverage  has  five  NH3  molecules  on  the 
surface  of  SWCNT  and  is  denoted  as  C5NH3.  Full  coverage  case  contains  ten  NH3 
molecules  and  is  expressed  as  C10NH3. 

The  structure  of  C5NH3  is  relaxed.  A  single  H2  is  attached  to  all  NH3  and  the 
structure  is  relaxed.  The  system  is  denoted  as  C5(NH3+H2).  The  binding  energy  is  0.10 
eV/H2.  Another  H2  is  attached  and  the  structure  is  relaxed.  The  binding  energy  per  H2  is 
0.097.  In  a  similar  way  we  have  attached  third,  fourth  and  fifth  molecule  on  each  NH3. 
All  the  constructed  systems  are  relaxed.  All  the  H2  are  physisorbed  and  are  equally 
perturbed  when  compared  to  the  unbound  H2  (H-H=0.74  A).  The  storage  capacity  of 
C5(NH3+5H2)  is  8.18  wt%  and  the  desorption  temperature  is  calculated  using  van’t  Hoff 
equation  and  is  1 15K. 

Naturally  we  can  expect  increase  in  storage  capacity  when  more  number  of 
adsorbates  is  introduced.  So  we  have  included  five  more  NH3  molecules  on  the  surface  of 
SWCNT.  The  clustering  of  NH3  molecules  is  noted.  On  addition  of  H2  the  clustering 
pattern  changed.  Initially  single  H2  is  attached  to  each  NH3  in  C10NH3.  The  structure  is 
relaxed  and  the  binding  energy  per  H2  is  0. 1  eV.  After  that  second,  third,  fourth  and  fifth 
H2  are  attached  and  the  structures  are  relaxed.  The  binding  energy  is  noted  in 
C10(NH3+4H2)and  is  0.098eV/H2.  The  corresponding  desorption  temperature  is  125K. 
The  maximum  storage  capacity  corresponds  to  C10(NH3+5H2)  and  is  13.2  wt%.  Since 
the  clustering  is  weak  it  has  no  effect  on  the  storage  capacity.  But  the  bonding  nature  of 
some  adsorbed  H2  is  changed. 

We  conclude  that  functionalization  with  electron  rich  and  neutral  molecule,  NH3 
enhances  the  storage  capacity  of  SWCNT  to  a  larger  extent.  The  uniform  perturbation  of 
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H-H  bond  lengths  in  the  case  of  half  coverage  is  noted  and  it  is  a  new  feature.  Also  we 
can  infer  that  both  in  half  coverage  and  full  coverage  cases  the  binding  energy  lie  below 
the  range  recommended  for  efficient  hydrogen  storage  medium.  So  due  to  the  low 
binding  energy,  H2  desorption  takes  place  below  the  room  temperature. 

(ii)  Aluminum  hydride  coated  single  walled  carbon  nanotube  as  a  hydrogen  storage 

medium 

We  have  taken  the  (5,  5)  single  walled  carbon  nanotube  for  our  study.  The  AIH3 
molecule  is  allowed  to  approach  the  CNT  from  large  distances,  to  determine  the 
energetically  preferred  absorption  position  of  AIH3.  The  AIH3  molecule  is  attached  to  the 
SWCNT  and  this  system  is  denoted  as  CAIH3.  The  binding  energy  of  AIH3  with  CNT  is 
found  to  be  1.46eV.  This  serves  as  the  reference  system  (CAIH3).  Interestingly,  AIH3 
group  forms  complex  with  hydrogen  molecules.  To  find  out  the  different  reaction 
products  between  H2  and  CAIH3,  a  series  of  energy  minimization  calculations  with  one, 
two,  three,  four,  five  and  six  H2  molecules  and  optimization  of  the  atomic  positions  (x,  y, 
z)  have  been  carried  out.  Hydrogen  molecule  is  allowed  to  move  freely  from  large 
distance  with  the  initial  bond  length  of  0.747A.  It  moves  towards  CAIH3  and  interacts 
with  the  system  and  gets  adsorbed  with  an  increased  bond  length  of  0.768A.  We  obtained 
the  binding  energy  as  0.16  eV  for  the  first  H2.  There  is  also  a  possibility  to  attach  more 
number  of  H2  to  this  system.  In  the  second  H2  adsorption,  the  H-H  distance  is  0.749A  and 
it  is  slightly  low  compared  to  the  first  H2.  The  binding  energy  is  reduced  to  0.09eV. 

So  there  is  a  possibility  to  attach  more  number  of  AIH3  to  the  nanotube  to  increase 
the  hydrogen  storage  capacity.  We  have  investigated  the  two  cases  namely  C5AIH3  and 
CIOAIH3. 

As  a  result  each  AIH3  can  bind  up  to  three  or  four  hydrogen  molecules.  These 
systems  are  quite  stable  and  the  H2  can  be  extracted  from  the  system  without  disturbing 
the  C-Al  bonding  or  detaching  the  A1H3  from  the  carbon  nanotube.  Of  the  two  systems 
C5AIH3  (half)  and  CIOAIH3  (full),  C5AIH3  system  is  qualified  as  a  viable  hydrogen 
storage  material.  The  average  binding  energy  per  H2  is  in  the  range  of  0.10  eV  to  0.2  eV 
which  is  recommended  for  a  HSM.  The  hydrogen  storage  capacity  for  half  coverage  is 
8.3  wt%.  This  value  meets  the  standard  fixed  by  U.S.  Department  of  Energy.  Here  again 
the  H2  release  is  below  the  room  temperature. 

(iii) Functionalization  of  single  walled  carbon  nanotube  with  borane  for  hydrogen 

storage 

We  have  chosen  another  light  weight  molecule  BH3.  It  is  also  hydrogen  rich.  But  it  is 
electron  deficient.  This  study  highlights  the  high  adsorption  capability  of  this  electron 
deficient  adsorbate  dispersed  on  the  surface  of  the  CNT.  The  half  coverage  case  (C5BH3) 
has  five  BH3  molecules  and  the  full  coverage  case  (CIOBH3)  has  ten  BH3  molecules  on 
the  full  hexagons  of  CNT. 

Initially  single  H2  is  attached  to  each  BH3  in  the  relaxed  structure  of  C5BH3.  This 
structure  is  relaxed  and  the  H2  has  been  displaced  close  to  BH3.  Similarly  second,  third, 
fourth  and  fifth  H2  molecules  are  attached  to  each  BH3.  The  systems  are  relaxed.  The 
bond  lengths  noted  after  relaxation  reveal  that  one  of  the  H2  is  bound  closely  to  BH3 


which  has  resulted  in  the  H-H  bond  lengthening.  This  feature  is  noted  in  all  the  systems. 
If  the  desorption  has  to  take  place  at  ambient  conditions,  the  ideal  H2  binding  energy 
range  should  be  0.2  -  0.4  eV.  The  system  C5(BH3+4H2)  has  the  binding  energy  0.21 
eV/H2.  This  lies  in  the  energy  range  required  for  a  viable  HSM.  The  storage  capacity  is 
6.83  wt%.  It  is  almost  equal  to  the  target  set  by  U.S.  DOE.  The  corresponding  desorption 
temperature  is  calculated  using  van’t  Hoff  equation  and  it  is  268K. 

We  intend  to  increase  the  storage  capacity  of  the  system  by  introducing  five  more 
BH3  molecules.  The  system  is  denoted  as  CIOBH3.  Similar  to  the  half  coverage  case  we 
have  attached  hydrogen  molecules  to  each  of  the  BH3  molecules  in  the  relaxed  structure 
of  CIOBH3  one  by  one,  up  to  a  maximum  of  five  H2.  The  clustering  started  in 
C10(BH3+3H2).  Two  of  the  BH3  molecules  have  come  closer  at  a  distance  of  2.13  A 
between  them.  The  binding  energy  per  H2  in  C10(BH3+3H2)  is  0.26eV.  The  desorption 
temperature  is  328K  and  is  well  above  the  room  temperature.  The  storage  capacity  is 
8.9wt%.  A  single  dimer  is  observed  in  C10(BH3+4H2)  and  C10(BH3+5H2).  Some  of  the 
H2  molecules  are  dissociated  into  H  atoms.  We  can  expect  these  H  atoms  to  combine  as 
H2  while  desorption.  Remaining  H2  are  molecularly  physisorbed.  The  binding  energy  per 
H2  in  C10(BH3+4H2)  is  0.24eV.  The  corresponding  desorption  temperature  is  306K  and 
is  nearly  equal  to  the  room  temperature.  The  storage  capacity  of  this  system  is  1 1.5wt%. 

We  conclude  that  light  weight,  hydrogen  rich  and  electron  deficient,  BH3  enhances 
the  hydrogen  storage  capacity  of  SWCNT  and  this  system  has  notable  binding  energy 
which  lies  in  the  range  recommended  for  achieving  DOE  target  at  ambient  conditions. 
The  systems  C5(BH3+4H2)  and  C10(BH3+4H2)  are  capable  of  releasing  H2  around  room 
temperature.  So  the  SWCNT  functionalized  with  BH3  becomes  as  a  practicable  hydrogen 
storage  medium  for  fuel  cells. 

(iv)  Molecular  hydrogen  storage  in  SWCNT  functionalized  with  Nickel  hydride 

In  this  study  (5,  5)  SWCNT  is  functionalized  with  Nickel  hydride  (NiH2)  molecules. 
The  simulation  cell  contains  two  unit  cells  of  (5,  5)  SWCNT.  The  half  coverage 
comprises  five  NiH2  molecules  attached  to  the  relaxed  CNT.  Initially  single  hydrogen 
molecule  is  attached  to  every  single  NiH2  and  after  relaxation  the  binding  energy  per  H2 
is  calculated  as  0.059eV.  On  subsequent  addition  of  H2  molecules  a  gradual  decrease  in 
the  binding  energy  per  H2  is  noted.  For  C5(NiH2+4H2)  an  increase  in  binding  energy  per 
H2  is  noted  and  is  0.1 4eV.  Some  of  the  attached  H2  molecules  are  absorbed  by  MH2 
molecules,  which  have  resulted  in  the  increase  of  the  binding  energy.  The  binding  energy 
per  H2  in  C5(NiH2+5H2)  is  0.079eV.  The  desorption  temperature  for  hydrogen  release  is 
calculated  as  175K  and  101K  for  C5(NiH2+4H2)  and  C5(NiH2+5H2).  The  storage 
capacities  of  C5(NiH2+4H2)  and  C5(NiH2+5H2)  are  4.89  wt%  and  6.05  wt%. 

The  full  coverage  encloses  ten  NiH2  molecules  in  (5,  5)  SWCNT  and  the  composition 
is  relaxed.  The  clustering  of  N1H2  molecules  is  observed  in  contrast  to  the  half  coverage 
where  the  occurrence  is  absent.  On  further  addition  of  H2  molecules  the  clustering  pattern 
is  found  to  change.  A  gradual  decrease  in  the  binding  energy  per  H2  is  noticed.  The 
system  C10(NiH2+H2)  and  C10(NiH2+2H2)  have  the  binding  energy  0.19eV/H2  and 
0.12eV/H2  ,  correspondingly  the  desorption  temperatures  are  246K  and  157K.  The 
binding  energy  for  single  H2  in  C10(NiH2+4H2)  and  C10(NiH2+5H2)  are  0.082eV  and 
0.056eV  respectively.  The  corresponding  desorption  temperatures  are  104K  and  72K. 


Due  to  clustering  the  hydrogen  atoms  in  NiH2  have  come  out  and  some  has  combined  as 
H2  molecules.  This  feature  is  detected  in  all  configurations  of  full  coverage.  We  can 
expect  these  molecules  obtainable  for  use  in  the  fuel  cells.  So  the  clustering  has  enhanced 
the  storage  capacity  of  the  system.  The  storage  capacities  of  C10(NiH2+4H2)  and 
C10(NiH2+5H2)  are  6.91  wt%  and  8.50  wt%  .  The  binding  energy  per  H2  in  all  the 
systems  indicates  that  all  H2  are  physisorbed.  Since  the  binding  energy  per  H2  gets 
reduced  on  higher  coverage  of  H2  molecules,  the  systems  deliver  hydrogen  below  the 
room  temperature. 

Summary: 

Out  of  the  four  systems  investigated,  the  single  walled  Carbon  Nanotube 
functionalized  with  BH3  is  found  to  be  suitable  hydrogen  storage  medium.  It  has 
hydrogen  storage  capacity,  1 1.5  wt%  and  the  corresponding  H2  desorption  temperature  is 
306  K.  The  system  has  to  be  fine  tuned  for  a  better  performance. 
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Work  in  Progress  with  the  follow  on  project: 

1.  Study  of  the  effect  of  finite  length  of  SWCNT  hydrogen  storage  using  metal 
hydrides 

2.  Molecular  dynamics  calculations  for  all  the  systems  investigated  are  to  be 
completed. 

3.  The  SWCNT  (5,  5)  is  functionalized  with  lithium  hydride,  LiH  and  its  storage 
capacity  is  being  investigated. 

4.  To  study  the  effect  of  curvature  of  the  tube  on  the  hydrogen  storage,  we  are 
investigating  large  diameter  tubes  (10,0)  and  (10,10). 

5.  Searching  for  other  materials  which  will  be  capable  of  releasing  H2  at  ambient 
conditions. 


Hydrogen  Storage  in  BH3  and  NH3 
functionalized  (5,5)  SWCNT 


Hydrogen  storage  in  CNTs: 


•  The  high  surface  to  volume  ratio  ,  light  weight  and 
extensive  gas  adsorption  properties  made  CNTs  to  be 
considered  as  the  promising  candidates  for  hydrogen 
storage. 

•  The  experimental  results  by  Dillon  et  al  formed  the  basis 
for  the  hydrogen  storage  in  CNTs.  Since  then  a  lot  of 
works  have  been  devoted  to  explore  the  hydrogen 
storage  capacity  of  CNTs. 

•  The  interaction  between  H2  and  bare  CNT  is  weak  and 
so  it  can't  remain  attached  even  at  ambient  conditions. 
So  functionalization/  doping  of  CNTs  with  various  atoms, 
molecules  and  metal  complexes  enhances  the  storage 
capacity  when  compared  to  bare  CNTs. 


The  functionalization  of  CNT  with  transition  metal 
atoms  has  increased  the  storage  capacity  to  a 
higher  extent  than  the  value  predicted  by  U.S. 
Department  of  energy.  But  transition  metal  have 
heavy  mass. 

To  achieve  high  gravimetric  and  volumetric 
densities  it  is  important  to  find  light  weight 
atoms/ complexes  which  enhances  the  storage 
capacity  of  SWCNT. 

We  have  chosen  BH3  and  NH3  molecules  for 
functionalization  and  they  have  molecular  masses, 
13.8  gmol 1  and  17.03  gmol 1  respectively 

Both  being  light  weight  molecules  and  hydrogen 
rich  can  adsorb  hydrogen  to  certain  extent  and  our 
study  highlights  this  fact. 


bh3 

•  Borane  or  Boron  trihydride  has  molecular  structure 
analogue  to  AIH3. 

•  It  exists  in  gaseous  form. 

•  It  is  the  parent  compound  of  all  boranes. 

•  Generally  boranes  are  electron  deficient  and  highly 
reactive. 


nh3 

•  Ammonia  or  Azane  exists  in  gaseous  form. 

•  1 1  has  trigonal  pyramidal  shape. 

•  It  has  lone  pair  of  electrons. 


Methodology  and  modelling 


The  total  energy  calculations  are  carried  out  using  density  functional 
theory  implemented  in  VASP  code. 

Two  unit  cells  of  SWCNT  (5,5)  in  the  simulation  cell  of  dimension 
24x24x4. 98A3  is  taken. 

Ultra  soft  pseudopotentials  are  used  with 

C:2s2  2p2 

B:  2s2  2px 

N:2s2  2p3 

H:ls 

as  valence  states. 

The  coordinates  of  SWCNT  (5,5)  are  generated.  BH3  and  NH3 
molecules  are  attached  to  the  relaxed  geometry  of  CNT,  resulting  in 
the  formation  of  CBH3  and  CNH3  systems. 

The  investigation  is  done  for  half  coverage  which  has  5  molecules  and 
full  coverage  which  has  10  molecules  of  BH3  and  NH3. 


C5BH3 

The  hydrogen  molecules  are  attached  to  the  relaxed 
structure  of  C5BH3.  The  structures  are  optimized.  The 
binding  energy  per  H2  and  the  storage  capacities  are 
calculated. 


System 

Binding  energy 
per  H2  (eV) 

Storage  capacity , 
wt% 

C5(BH3+H2) 

0.65 

1.8 

C5(BH3+2H2) 

0.35 

3.53 

C5(BH3+3H2) 

0.26 

5.21 

C5(BH3+4H2) 

0.21 

6.83 

C5(BH3+5H2) 

0.18 

8.4 

C5(BH3+H2) 


C5(BH3+3H2) 


C5(BH3+4H2) 


C10BH3 

As  in  the  case  of  half  coverage  five  hydrogen  molecules 
are  added  one  by  one  to  each  BH3  in  the  relaxed 
structure  of  C10BH3.  The  binding  energy  for  single  H2 
and  storage  capacity  of  all  relaxed  systems  are 
calculated. 


System 

Binding  energy  per 
H2,eV 

Storage 
capacity,  wt% 

C10(BH3+H2) 

0.67 

3.15 

C10(BH3+2H2) 

0.36 

6.12 

C10(BH3+3H2) 

0.26 

8.90 

C10(BH3+4H2) 

0.24 

11.5 

C10(BH3+5H2) 

0.18 

14 

C10(BHo+H2) 


C10(BH3+3H2) 


C10(BH3+4H2) 


C5NH3 

As  in  C5BH3,  five  hydrogen  molecules  are  attached  to 
each  NH3  in  C5NH3  .  The  binding  energy  per  H2  and 
storage  capacity  of  all  optimized  systems  are 
calculated. 


System 

Binding  energy  per 
H2,eV 

Storage 
capacity  wt% 

C5(NH3+H2) 

0.10 

1.75 

C5(NH3+2H2) 

0.097 

3.44 

C5(NH3+3H2) 

0.095 

5.07 

C5(NH3+4H2) 

0.092 

6.65 

C5(NH3+5H2) 

0.091 

8.18 

C5NH3 
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C5(NH3+H2) 


<5> 

Q)  9) 


8 

T  Q) 
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<9 

Q) 


C5(NH3+4H2) 


C5(NH3+5H2) 


C10NH3 

On  addition  of  another  H2  to  C10(NH3+h2)  ,  we 
observed  decrease  in  the  binding  energy  and  it  is  low. 
The  binding  energy  per  H2  and  storage  capacities  of  all 
the  systems  are  calculated. 


System 

Binding  energy  per 
H2,eV 

Storage 
capacity  wt% 

C10(NH3+H2) 

0.10 

3 

C10(NH3+2H2) 

0.099 

5.83 

C10(NH3+3H2) 

0.097 

8.50 

C10(NH3+4H2) 

0.098 

11.1 

C10(NH3+5H2) 

0.092 

13.2 

The  ideal  H2  binding  energy  range  to 
achieve  DOE  target  is  0.2-0.4eV  and  so 
desorption  takes  place  at  ambient 
conditions. 

Our  study  reveals  that  C5NH3+nH2  has  less 
binding  energy  when  compared  to  that  of 
C5BH3+nH2  (n=l,2,3,4&5). 


cbh3 


C5BH3 

•  On  addition  of  hydrogen  molecules  one  of  the  B-H  bond  length 
has  increased  from  1.18A  to  1.20A.Since  BH3  is  electron 
deficient  there  is  a  significant  charge  transfer  between  BH3  and 
H2  molecules. 

•  Moreover  in  each  system  one  hydrogen  molecule  is  closely 
bound  to  BH3  at  a  distance  of  1.32A.  The  bond  length  has 
extended  from  0.74A  to  0.85A. 

•  Remaining  H2  molecules  are  less  perturbed  and  their  bond 
length  has  increased  to  0.749A. 

C10BH3 

•  Here  too  in  full  coverage  one  of  the  adsorbed  H2  molecules  is 
more  perturbed  and  the  bond  length  is  0.85A  . 

•  Single  BH3  dinner  is  observed  in  C10(BH3+4H2)  and 
C10(BH3+5H2).The  formation  of  the  dimer  doesn't  seem  to 
enhance  the  binding  energy  per  H2  molecule. 


cnh3 


C5NH3 

•  All  the  hydrogen  molecules  perturbed  equally.  The  H-H 
bond  distance  has  increased  uniformly  to  0.75A. 

•  A  gradual  decrease  in  the  binding  energy  per  H2  is 
noted  when  compared  to  C5BH3+nH2 

( n=l,  2, 3, 4&5)  systems. 

CIONH3 

•  As  in  C5NH3  the  hydrogen  molecules  are  almost 
equally  perturbed. 

•  Here  clustering  between  NH3  molecules  is  noted. 

Since  the  clustering  weak  it  doesn't  increase  the 
binding  energy  of  single  H2. 


Salient  Features: CBH3 

The  system  C5(BH3+4H2)  has  the  binding 
energy  per  H2  as  0.21eV.The  corresponding 
desorption  temperature  is  calculated  using  van't 
hoff  equation  and  it  is  268K. 

The  system  C10(BH3+4H2)  has  a  notable  binding 
energy  of  0.24eV.The  desorption  temperature  is 
306K  and  it  is  equal  to  room  temperature.  The 
storage  capacity  exceeds  the  DOE  target  and  is 
11.5wt%. 


Salient  features  :  CNH3 

The  binding  energy  per  H2  in 
C5(NH3+5H2)  is  0. 092eV.The  desorption 
temperature  is  116K.  1 1  is  far  below  the 
room  temperature. 

In  C10(NH3+2H2)  the  binding  energy  for 
single  H2  is  0.099eV  and  the  desorption 
temperature  is  estimated  as  126K. 


Conclusion 


•  On  comparing  both  the  systems,  CBH3  qualifies 
as  an  viable  HSM  with  high  hydrogen  storage 
capacity  and  desorption  of  hydrogen  takes  place 
at  room  temperature. 

•  The  functionalization  with  electron  deficient  light 
weight  molecules  like  BH3  enhances  the 
hydrogen  storage  capacity  to  a  larger  extent. 

•  The  functionalization  with  electron  rich  molecule 
like  NH3  enhances  the  storage  capacity  but  the 
binding  energy  per  H2  gets  decreased  and  so 
they  can't  release  H2  at  room  temperature. 
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